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ABSTRACT 



Although only 22 millisecond pulsars (MSPs) are currently known to exist in 
the globular cluster (GC) 47 Tucanae, this cluster may harbor 30 — 60 MSPs, or 
even up to ~ 200. In this Letter, we model the pulsed curvature radiation (CR) 
gamma-ray flux expected from a population of MSPs in 47 Tucanae. These MSPs 
produce gamma-rays in their magnetospheres via accelerated electron primaries 
which are moving along curved magnetic fleld lines. A GC like 47 Tucanae con- 
taining a large number of MSPs provides the opportunity to study a randomized 
set of pulsar geometries. Geometry-averaged spectra make the testing of the 
underlying pulsar model more reliable, since in this case the relative flux uncer- 
tainty is reduced by one order of magnitude relative to the variation expected 
for individual pulsars (if the number of visible pulsars = 100). Our predicted 
spectra violate the EGRET upper limit at 1 GeV, constraining the product of 
the number of visible pulsars A^ and the average integral flux above 1 GeV per 
pulsar. GLAST/LAT should place even more stringent constraints on this prod- 
uct, and may also limit the maximum average accelerating potential by probing 
the CR spectral tail. For A^ = 22 - 200, a GLAST/LAT non-detection will lead 
to the constraints that the average integral flux per pulsar should be lower by 
factors 0.03 — 0.003 than current model predictions. 

Subject headings: gamma rays: theory — globular clusters: individual (47 Tu- 
canae) — pulsars: general — radiation mechanisms: non-thermal 
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Introduction 



Millisecond pulsars (MSPs) h ave been detected in several globula r clusters (GCs), 
including the nearby 47 Tucanae ( Camilo et al. 200ol : Freire et al. 2003 ). Although only 



22 M SPs are currently known to exist in 47 Tucanae (for a recent re yiew. see iLorimer et al. 



20031). it is expect ed that this cluster m ay harbor 30 — 60 MSPs (ICarailo fc Rasio 
Heinke et al. 2005 ). or even up to ~ 200 ( Ivanova. Fregeau. fc Rasio 20051 ). 



2005 



In this Letter, we model the pulsed curvature radiation (CR) gamma-ray flux expected 
from a population of MSPs in 47 Tucanae (§ [2]), and compare it with the GLAST/LAT 
sensitivity and EGRET upper limits. Our approach differs from similar recent s tudies 
JSednarek k Sitarekl I2OO7I : iHarding. Usov. k Mushmovl I2OO5I : Icheng k Taaml I2OO3I . here- 
after BS07, HUM05, k CT03) in a number of respects. We firstly use a fully 3 D General- 
Relativistic (GR) polar cap (PC) pulsar model (see e.g. I Venter k de Jagerll2005l ). sampling 
non-thermal gamma-ray radiation from a range of magnetic field lines above the PC. Next, 
our study involves a population of MSPs, and we sample randomly from this ensemble. 
We are therefore able to calculate absolute cumulative fiuxes expected from the MSPs in 
47 Tucanae, without needing to use approximate spectra, or scale up single MSP predictions 
or single particle spectra. We are also specifically interested in the range of pulsed spec- 
tra expected due to the uncertain geometries and spread in values of the period P and its 
derivative P of these pulsars. 

X-rays from the MSPs in 47 Tucanae are believed to be of thermal orig i n, and are likely 

prod uced by a return current heating the pulsars' PCs f Bogdanov et al.|[2006 : Harding k Muslimov 
2OO2I . CT03). These X-rays may be Compton upscattered to high-energy gamma-rays. How- 
ever, this inverse Compton scatterin g (ICS) component is believ ed to be dominated by the 
CR component (see e.g. Figure 3 of iBulik. Rudak. k Dyksll2000l ). and we therefore omitted 
it from our calculations below. 

We lastly make some conclusions regarding the possibility of constraining a model of 
gamma-ray radiation from GCs (§11]), since an ensemble of MSPs in a GC cluster provides 
a unique opportunity for sp ecifically constraining the GR frame-dragging pulsar model (e.g. 
Harding k Muslimovlll998l ) in a geometry- independent way (§[3]). 



2. Pulsed Gamma-Ray Flux 



The measured values of P for GC pulsars are affected by their acceleration in the 
gravitational potential of the GC, and therefore need to be corrected to ob tain intrinsic 
period derivatives Pint- This was done for 47 Tucanae by Bogdanov et al. (j2006 ) who derived 
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spin-down luminosities E^-ot assuming a King model. From their Table 4, we selected 1 3 MSP s 
and calculated their corresponding Pint using their periods P as given by lFreire et al.l (120031 ) . 
and the expression i^rot = 47r^/Ns-Pint/-P^5 with I^s the moment of inertia. 

The details of the implementatio n of our isolated MSP mod el, using the GR- framework 

of Har ding, Muslimov, & Tsygan (e.g. lHarding fc Muslimovll 19981 ) may be found in I Venter fc de Jager 
( l2005l ). where we discuss the essential link between pulsar visibility and the assumed geom- 
etry, i.e. magnetic inclination angle x observer angle ( (both measured with respect to 
the spin axis fi). 

HUM05 found that most MSPs are inefficiently screened by CR and ICS pairs. Screen- 
ing caused by electron-positron pairs will lower the accelerating potential, leading to a 
lower pulsed CR spectral cut-off. Higher numbers of energetic electrons may however 
boost the unpulsed gamma-ray spectrum. The vast majority of MSPs in our population 
l ie below the critical spin-down lu minosity above which screening is expected to occur 
faarding. Muslimov. fc Zhanejl2002[ l 



rot, break 



1.4 X 10 



34 



P 

BI2 



1/7 



ergs/s. 



(1) 



Here, B12 = G, and P is the pulsar period in seconds. (Note the positive sign of the 

power of P). This condition may be rewritten so that screening occurs when 



logio P > 2-625 logio P - 12.934 



(2) 



assuming /^s = OAMR"^, stellar radius R = 10^ cm, and pulsar mass M = 1.4M0. We used 
an unscreened electrical potential fo r the bulk of the popul ation, and used the approximation 
of a screened electric field given by iDyks &: RudakI (I2OOOI ) for MSPs with spin-down powers 
greater than this critical spin-down power (which was only the case for 47 Tuc U)). 

For each pulsar i in our populatio n, we calculated phase-averaged CR spectra via (see 
also Eq. 12 of 



Venter fc de Jagerll2005l ) 



dE 



1 



X 



sin (d(dE^ 




^, E^ + dE^) 



ES 



dL\{(t)L,xX,E^ 
d(l)Ld(dEy 



(3) 



bidCdE^, (4) 



where (pi is the phase angle, and I{E^, E^ + dE^) picks out the gamma-ray energies in the 
range {E^,E^ + dE^). Primary electrons leave a stellar surface patch dS relativistically 
(/3o = Vq/c ~ 1) at a rate of dN^ = pedSPoc/e, where pe is the charge density. These 
primaries radiate CR at a position (r, 6, 0) above the PC with a characteristic energy of 
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E!^^ = 1.5\cj^meC^p~^ , where Ac = h/{mec) is the Compton wavelength, and pc{r,9,(f)) is 
the curvature radius of the GR-corrected dipolar magnetic field. The CR is associated 
with an incremental gamma-ray luminosity of dL^ = dN^EQ-^dt, which is radiated in a 
time dt. The CR loss rate is given by Eck ~ 2e^c7^/(3p^), where 'y{r,9,(j)) is the Lorentz 
factor of the electron primaries. We calculated CR photon spectra dN^/dE{E^,x,C) for 
X = C = 10°, 20°, ■■■,80°, and i = 1,2, ...,13, assuming R = 10^ cm, M = IAMq, and 
/ns = 0.4M-R^ All spectra (13 x 8 x 8 = 832) were scaled to a distance of c? = 5 kpc 
(ICratton et"aDl2003h . 

We next randomly chose = 100 visible MSPs (with random x and (), and summed 
their pulsed spectra to obtain a million cumulative spectra from the Monte Carlo process: 

/ cum k=l 

where for each index j = 1,2, ...,Nt = 10^, a total of = 100 spectra (randomly sampled 
from the abovementioned 832 spectra) were summed. Therefore, j corresponds to a specific 
choice of x, C, and fc- values when oversampling from 13 to 100 pulsar spectra. This procedure 
was repeated for 1 MSP instead of 100, i.e. setting A^ = 1. Upon comparison of these results, 
one would expect that the mean values of the single MSP differential spectra would scale 
with A^, and the standard deviations with \/N. In panel (a) of Figure [H the mean values 
of E'^dN^I dE at three different photon energies E^ are shown as a function of the number 
Nt of values used to calculate the mean (with a maximum of A"* = 10^). The thick lines 
indicate results when A^ = 100, while the thin lines indicate results for A^ = 1, scaled by a 
factor 100. Similarly, panel (b) of Figure [T] indicates the standard deviation of E^dN^/dE at 
three different photon energies vs. A^t. The thick lines again indicate results when A^ = 100, 
while the thin lines indicate results for A^ = 1 scaled by a factor 10. It is clear that the mean 
and (T-values converge beyond A^ ~ 10^, for both A^ = 1 and A^ = 100. The A^ = 1 results 
are more erratic at small Aj-values. The relative error on the mean values are therefore a 
factor 10 lower when using a population of 100 MSPs, as opposed to the much larger error 
obtained for single pulsars. Figure [2] depicts average cumulative pulsed differential spectra 
(E^dN^/dE)^^^ for A^ = 100 and Nt = 10^, as well as 2o--bands. Also shown are the 
GLAST/LAT "sensitivity curve (HUM05) and EGRET upper hmits Jpierro et al.lll995h . as 
well as other predictions of pulsed gamma-ray radiation from 47 Tucanae (HUM05). 



3. Discussion 



The gamma-ray visibility of an isolated MSP depends extremely sensitively on the as- 
sumed pulsar geometry. Inclination and observer angles x and ( are usually estimated 
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from radio polarization measurements (e.g. [Manchester fc Johnston! (jl995l )). involving the 



rotating vector model. In our MSP model, the el ectric potential boundary condition of 



= 1) = is used (IMuslimov fc Harding 119971 ) as the last open magnetic field lines 



= 1) are treated as equipotentials. In addition, the space-charge-limited nature of this 
model requires that the accelerating electric field parallel to the magnetic field {E\\) should 
be zero at the stellar surface {r = R). This means that on-beam radiation (when x ~ C) will 
have a spectral cut-off which roughly scales with the maximum potential $max- Off-beam 
radiation will however rapidly decrease as the impact angle P = C ~ X increases, because 
this radiation is produced along field lines where the electric potential $, and hence E\\, has 
dropped significantly. 

An ensemble of MSPs provides the unique opportunity to side-step the issue of pul- 
sar geometry (which is crucial when modeling a single MSP) by averaging over numerous 
spectra and geometries, which results in a "geometry-averaged" spectrum with relatively 
small uncertainty. We also exploit this scenario in the sense that we use results from the 
corotating frames of the MSPs, since we expect that the effect of Lorentz transformations 
to the observer frame will average out over the ensemble. We find that the relative spread of 
ensemble values for the bolometric gamma-ray and electron luminosities as well as the spin- 
down power and average gamma-ray spectra are typically one order smaller than for the same 
quantities averaged as single MSP quantities (i.e. they scale with N~^/'^). This illustrates 
the importance of testing pulsar models by observing GCs. (This conclusion follows when 
using a constant equation of state (EOS). Variations in the EOS have a much smaller effect 
than that of differing pulsar geometries, and is therefore not included in this discussion; see 
§4). 

Due to the lower relative ensemble errors, one should therefore expect to be able to 
constrain average pulsar parameters using average spectra and their errors shown in Figure [2J 
This would not be the case if ensemble values did not converge (Figure [T]) and had relatively 
small errors. Furthermore, the fact that the ratios of average cumulative luminosities and 
spin-down powers to their single MSP counterparts are very close to = 100 gives us 
confidence in our sampling algorithm. 

In our model, the predicted average single pulsar efficiencies for electron and gamma-ray 
production are typically ~ 2% and ~ 7%, depending on the population of MSPs used. For 
our 47 Tucanae MSP population (with (^Erotj = 2.2 x 10^^ ergs/s), we found 0.74% and 
6.1% respectively. This should be compared to the esti mate of BS07 of 1% for the first 



value, and to that of iHarding. Muslimov. fc Zhand (120021 ) of ~ 10% for the latter quantity. 



Venter fc de Jagerl (120051 ) quoted values of 1 — 2.5% and 2 — 9% for these quantities for the 



case of the pulsar PSR J0437-4715. 
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In Figure [21 we converted integral EGRET upper limits (jPierro et al.lll995l ) and model 
predictions of HUM05 to differential values assuming a E'"^ spectrum. Our predicted pulsed 
gamma-ray flux at 100 MeV is below the differential EGRET upper limit, much lower than 
that of HUM05. Their prediction of the integral flux above 100 MeV exceeds the EGRET 
upper limit by a factor / ~ 19 for 15 MSPs, and / ~ 125 for 100 MSPs. (T his seemingly over- 
predi ction of integral flux was also found in the case of PSR J0437-4715 (IVenter fc de Jager 
20051 ). HUM05's prediction is a factor / ~ 35 above the EGRET upper limit for this 
MSP). However, o ur summed pulsed spectra violate the EGRET upper limit at 1 GeV 
(IFierro et al.lll995l ). providing constraints on the low-energy tail of the average pulsed spec- 
trum. GLAST/LAT observations will however constrain the predicted spectrum at all en- 
ergies of relevance, which should provide a meaningful constraint on the product of visible 
pulsars N, and the average intergral flux above 1 GeV per pulsar. 



4. Conclusions 

In this Letter, we modeled the pulsed gamma-ray flux expected from a population of 
MSPs in 47 Tucanae. We chose a number of = 100 visible members, but the spectra can 
easily be (linearly) scaled to any other reasonable number. We were especially interested in 
obtaining errors on the pulsed gamma-ray spectrum which would reflect the uncertainty in 
the geometries, as well as the spread in P and P, of the individual pulsars. In this way, we 
wanted to see if it would be possible to constrain the underlying pulsar model independently 
from the assumed pulsar geometry. It remains to be shown if the spread of the 13 selected 
pulsars represent the true inherent spread of MSPs in 47 Tucanae. 

For completeness' sake, we furthermore investigated the effect of varying the EOS (M, R, 
and /). By choosing Jns = OAMR'^, and letting M and R vary between M/Mq = 1.4, 1.5, 1.6 
and Rq = R/10^ cm = 1.0, 1.2, 1.4, 1.6, we found that the maximum integral flux did not vary 
by more than ~ 4% (although it increases with both M and R; the spectral cut-off energy 
furthermore did not vary by more than ~ 13%). Since this variation will not significantly 
influence our constraints derived below, we used fixed values of M = 1.4M0 and Rq = 1 
throughout. 

Our average pulsed spectrum, including errors, exceed the EGRET upper limit at 1 GeV 
(for the average integral flux spectrum, by a factor of ~ 3.9). Assuming that this model 
predicts the gamma-ray flux correctly, this EGRET upper limit constrains the number of GC 
pulsars to 100/3.9 ~ 25. Generally, the constraints derived may be summarised as follows: 

^ < 'A (6) 

\-r model/ Q 
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with (Fobs) and (Fmodei) the average observed and predicted integral flux above 1 GeV 
per pulsar, and g is a factor indicating telescope sensitivity. For EGRET, q = 1, while 
g = 40 for GLAST/LAT. E.g., setting N = 22, 100, and 200, we find that (Fobs) / (F^odei) < 
1.1,0.25,0.13 for EGRET. This indicates that the EGRET hmit is obeyed for = 22, 
but that larger values of A^ imply a reduction in the model-predicted integral flux (by e.g. 
factors of 4 and 8 for A^ = 100 and 200). In the case of a GLAST/LAT non-detection, we find 
(Fohs) I (i^modei) < 0.03, 0.006, 0.003 for A^ = 22, 100, and 200. GLAST/LAT might lastly 
even constrain the maximum average electric potential per pulsar ($max) via the CR cut-off 
energy. It should also in principle be possible to look for the radio periods corresponding to 
the 22 known MSPs in 47 Tucanae in the gamma-ray data (or detect new gamma-ray periods 
if beaming properties are different), should GLAST/LAT detect the cumulative pulsed CR 
spectrum. 

This publication is based upon research supported by the South African National Re- 
search Foundation and the SA Gentre for High Performance Gomputing. 
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Fig. 1. — The mean and standard deviation a of E^dN^/dE (in units of ergs/s/cm^) vs. A^t 
as shown in panel (a) and panel (b), with Nt the number of values of E^dN^/dE used to 
calculate the mean and cr-values. Line types correspond to different energies as shown in 
the legend. Thick lines are for values of E'^dN^/dE summed over N = 100 randomly chosen 
visible pulsars, while thin lines are for values of E^dN^/dE for single (iV = 1) randomly 
chosen pulsars. The latter values were scaled by a factor 100 in panel (a), and a factor 
= 10 in panel (b). Results converge beyond Nt ~ 10 . 
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Fig. 2. — Average differential pulsed gamma-ray spectrum (solid line) with 2cr-error bands 
(shaded band between dashed lines) for = 100 and Nt = 10^, with smoothed high-energy 
tails. Furthermore we indicate predictions of HUM05 (lower square for 15 MSPs, upper 
square scaled to 100 MSPs). The predictions from HUM05 in volves estimates of single- 
particle spectra scaled using the Goldreich- Julian PC current (iGoldreich &: JulianI Il969l ) 
and arbitrary beaming angles. Our spectra include variations across the PC (sampling 
different magnetic field lines) as well as variations due to various inclination angles and ob- 
server geometries (with beaming angles calculated im plicitly). The GLA ST/LAT sensitivity 
(HUM05) as well as EGRET upper limits (diamonds - Fierro et al. IQQsl ) are also shown (the 
EGRET upper limits and HUM05 results were converted to differential values assuming an 
£'~^-spectrum). 



